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Synergism  of  Analytical  and  Numerical 
Approaches  in  Combustion 

B.J.  Matkowsky 

Department  of  Engineering  Sciences 
and  Applied  Mathematics 
Evanston,  Illinois  60208,  USA 


We  employ  a  synergism  between  analysis  and  large  scale 
scientific  computations  to  investigate  problems  of  both 
solid  and  gaseous  fuel  combustion.  Since  activation  ener¬ 
gies  in  combustion  problems  are  generally  large,  reaction 
is  typically  confined  to  a  thin  zone.  In  the  asymptotic 
limit  of  very  large  activation  energies,  tho  reaction  zone 
shrinks  to  a  moving  surface,  termed  a  front,  and  the 
problem  is  reduced  to  a  free  boundary  problem,  in  which 
determination  of  the  front  is  part  of  the  solution  of  the 
problem.  Our  analytical  studies  correspond  to  the  case  of 
a  front,  while  our  numerical  computations  correspond  to 
finite,  though  large,  activation  energies.  Since  the 
solution  of  the  problem  exhibits  rapid  change  through  the 
reaction  zone,  accurate  resolution  of  solution  behavior 
poses  a  challenge,  which  is  met  by  an  adaptive  pseudo- 
spectral  method  which  we  introduce.  Our  investigations  are 
characterized  by  the  interplay  between  analysis  and  numeri¬ 
cal  computation. 

In  solid  fuel  combustion  a  high  temperature  thennal  wave 
propagates  through  a  sample,  converting  reactants  to 
products.  In  addition  to  the  uniform  mode  of  propagation, 
in  which  a  planar  front  moves  with  constant  velocity,  other 
modes  of  propagation  have  been  observed.  These  include  (i) 
pulsating  combustion  -  in  which  the  velocity  of  a  planar 
front  oscillates  in  time,  (ii)  spinning  combustion  -  in 
which  a  hot  spot(s)  moves  in  a  helical  path  along  the 
surface  of  a  cylindrical  sample,  (iii)  multiple  point 
combustion  -  in  which  the  hot  spots  appear,  disappear  and 
then  reappear,  and  (iv)  spiral  combustion  -  in  which  a  hot 
spot  moves  along  a  spiral  path  on  a  planar  sample.  We 
describe  these  and  other  combustion  modes,  including  those 
exhibiting  chaotic  behavior. 

In  gaseous  fuel  combustion  in  two  dimensions,  we  consider 
the  problem  of  a  flame  stabilized  on  a  fuel  source.  In 
addition  to  an  axisymmetric  flame  front,  separating  burned 
from  unburned  gases,  we  describe  various  types  of  flames, 
including  (i)  pulsating  axisymmetric  flames,  (ii)  statio¬ 
nary  cellular  flames,  (iii)  cellular  flames  with  traveling 
waves,  standing  waves,  guasi-periodic  waves  and  chaotic 
waves  along  the  front. 

Finally,  we  describe  a  method  which  combines  our  adaptive 
pseudo-spectral  method  with  domain  decomposition  ideas,  to 
solve  problems  with  more  complex  chemical  kinetics.  Speci¬ 
fically,  we  illustrate  the  method  for  a  problem  with  a 
sequential  reaction  mechanism. 
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Blow-up  of  Phenomena  for  Some  Parabolic  and  Hyperbolic 
Models  in  Combustion  Theory 


Alberto  Tesei 

Dipart imento  di  Matematica  "G.  Castelnuovo" 
Universita  di  Roma  "La  Sapienza",  Roma,  Italy 


Even  the  simplest  mathematical  models  of  combustion  pheno¬ 
mena  are  inherently  nonlinear.  As  a  consequence,  the 
superposition  principle  is  not  valid  and  novel  "dissipa¬ 
tive"  phenomena  arise,  due  to  the  mutual  interaction 
between  diffusion,  convection  and  bulk  heat  sources  or 
sinks . 

The  space-time  evolution  of  such  phenomena  can  be  qualita¬ 
tively  studied  in  detail  using  simple  parabolic  or  hyper¬ 
bolic  models.  A  number  of  methods  has  been  devised  for  such 
investigation.  We  mention  in  particular  the  group  analysis 
of  nonlinear  equations,  comparison  methods  and  the  con¬ 
struction  of  stable  numerical  algorithms  for  conservation 
laws  - 

In  this  lecture  we  review  some  of  the  above  topics,  with  an 
emphasis  on  phenomena  of  thermal  runaway.  In  this  general 
framework,  we  shall  also  present  new  results  concerning 
blow-up  of  solutions  for  quasilinear  hyperbolic  conserva¬ 
tion  laws. 


On  the  Role  of  Computer-Aided  Calculations 
in  Combustion  Theory 


A.G.  Merzhanov 

Institute  of  Structural  Macrokinetics 
Russian  Academy  of  Sciences 
142432  Chernogoloka ,  Moscow  Region,  Russia 


The  computer-aided  numerical  solutions  to  the  equations  of 
thermal  conductivity  and  diffusion  with  the  nonlinear  heat 
and  mass  sources  originating  from  chemical  reactions  were 
first  used  for  obtaining  the  quantitative  results  of  some 
specified,  particular  cases  of  interest.  Then  it  was 
understood  that  a  multitude  of  such  calculations,  termed 
the  "numerical  experiment",  could  provide  a  powerful  tool 
for  theoretical  studies. 

Discussed  is  the  role  of  computer-aided  calculations  in 
combustion  theory  aiming  at: 

-  Obtaining  the  interpolation  formulae  that  may  substitute 
the  explicit  solutions. 

-  Numerical  checking  the  approximate  explicit  solutions 
(provided  by  approximate  theories) . 

-  The  studies  on  the  mechanism  of  the  processes;  elucida¬ 
ting  the  major  and  minor  factors  (terms  of  equations) ; 
preparing  the  equations  for  finding  the  approximate 
solutions . 

Considered  is  a  number  of  theoretical  problems  solved 
by  the  author  and  his  colleagues  at  one  time  by  using 
computers  (hot  spot  thermal  explosion,  transition  from 
self-inflammation  to  ignition,  combustion  in  a  flow, 
degeneration  of  combustion  wave,  self-ignition  of  metal 
particles  and  suspended  drops,  auto-oscillating  combustion 
and  spinning  waves,  etc.). 


NONLINLAK  STABlLi'lY  JN  SOLID  PKOl’LLLAN'l'  LOMlu::;'!  1  ON : 

A  THLORLTJCAL  APPROACH 

C.D.Pagani,  S. Salsa,  M.V'erri 

Uiparl imontd  di  Ma I ema t i ca  del  Pol i tec n i ce ,  Milano,  Italv 


A  nonlinear  reaction-diffusion  equation  on  a  half-line,  arising  in 
solid  propellant  combustion  theory,  is  considered.  The  propellant  is  sub¬ 
jected  to  an  external  radiant  flux;  furthermore,  it  is  assumed  adiabatic 
and  optically  opaque  except  at  the  burning  surface,  where  chemical  reac¬ 
tions  arc  concentrated.  We  investigate  the  existence  of  uniformly  propa¬ 
gating  thermal  waves  and  discuss  their  asymptotic  stability.  The  problem 
is  formulated  in  an  abstract  setting  as  an  initial  value  problem  in  a 
suitable  Banach  space  for  a  fully  nonlinear  parabolic  equation.  The 
obtained  stability/instability  results  are  valid  for  a  large  class  of 
pyrolvsis  laws  and  flame  submodels  which  are  incorporated  in  the  governing 
equation  of  the  system. 
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<;n  theory  of  surface  spin  burning 


BORIS  V.NOVOZHILOV 

Institute  of  Chemical  Physics,  Russian  Academy  of  Sciences. 

Moscow,  117977,  Russia 


The  problem  of  spin  combustion  is  taken  as  an  example  to 
discuss  various  methods  (analytical,  numerical  and  approximate) 
in  combustion  theory. 

For  the  first  time  spin  burning  has  been  observed  during  the 
burning  of  porous  metallic  specimens  in  gaseous  oxidizer.  At 
ignition  of  the  cylinder  sample  a  glowing  spot  of  small  diameter 
appears  on  its  side  surface,  moves  in  a  spiral  along  the  side 
surface  of  an  unburnt  part  of  the  sample.  The  sample  burns  along 
the  side  svirface;  the  inner  part  of  the  sample  either  does  not 
burn  up  at  all,  or  reacts  completely  after  the  surface  combustion 
front  has  passed. 

The  spin  burning  problem  is  liable  to  analytical  solution 
only  in  the  vicinity  of  the  stability  boundary ,whe''e  the 
departure  of  the  burning  regime  from  the  steady-state  is 
considered  to  be  small.  A  number  of  papers  were  devoted  to  the 
theoretical  investigation  of  the  problem  with  the  help  of 
bifurcation  theory.  There  are  only  a  few  study  where  a  pronounced 
non-linear  regime  of  spin  combustion  was  treated  numerically. 

It  must  be  mentioned  that  the  spin  burning  has  been  realized 
both  analytically  and  numerically  with  very  small  cylinder 
diameters.  Real  experiments  are  run  ,as  a  rule,  with  specimens  of 
a  much  larger  diameters. 

The  approximate  analytical  method  is  formulated  for  the 
surface  spin  combustion  in  this  paper.  It  is  assumed  in  a 
two-dimensional  model  that  the  unburned  fuel  as  well  as  the 
combustion  products  are  both  solid  and  a  condensed  phase 
exothermic  reaction  occurs  on  the  cylindrical  surface,  with  the 
combustion  centre  penetrating  into  the  specimen  to  a  small  depth. 
Because  the  spin  spot  motion  is  almost,  perpendicular  to  the 
vector  of  the  combustion  wave  mean  velocity  the  two-dimensional 
problem  can  be  reduced  to  a  simple  one-dimensional  one,  which  is 
similar  to  the  problem  of  combustion  with  heat  losses. 

The  conditions  of  existence  of  spin  regime  are  obtained. 
Explicit  analytical  e-xpression  for  various  quantities  of 
intere.st,  including  the  propagation  spin  velocity,  the  wave  mean 
velocity  and  the  temperature  profile  are  given.  They  depend  on 
the  cylinder  dicvmeter  and  heat  losses.  As  the  diameter  grow.s  two 
or  more  centre.s  of  the  chemical  reaction  appear  ( multi -cent.re 
spin)  on  t}io  cylinder  surface. 

Wlierevei'  po.v.sible.  conclusions  of  the  above  mentioned 
approache.s  are  c;ompared  with  the  expei'i menta  1  data  and  the 
nunie  r  i  c  i\  1  re  su  1  tr. . 


On  a  Moving  Boundary  Problem  Arising 
in  Fluidized  Bed  Combustion 


R.  Dal  Passo 

Istituto  per  le  Applicazioni  del  CaJcolo  "M.  Picone" 
Consiglio  Nazionale  delle  Ricerche  Roma,  Italy 


We  consider  a  model  describing  the  combustion  of  a  coal 
particle  in  a  fluidized  bed,  in  which  attrition  plays  a 
dominant  role.  The  model  consists  of:  1)  a  quasilinear 
elliptic  equation  for  the  oxygen  concentration,  supple¬ 
mented  by  boundary  conditions  on  the  moving  surface  repre¬ 
senting  the  burning  particle  interface;  2)  an  ordinary 
evolution  equation  for  the  carbon  consumption;  3)  an 
equation  governing  the  motion  of  the  interface  in  terms  of 
a  specified  function  of  the  carbon  consumption  at  the 
interface.  We  prove  a  global  existence  and  uniqueness 
result,  together  with  a-priori  bounds  for  the  so..ution;  the 
existence  of  travelling  waves  will  also  be  established. 


NONLINEAR  TRANSIENT  BURNING  OF  SOLID  PROPELLANTS; 
A  SIMPLE  EXAMPLE  OF  SYNERGISM 
BETWEEN  ANALYTICAL  AND  NUMERICAL  METHODS 


Luigi  Dc  Luca 

Dipartimento  di  Energetica,  Politecmco  di  Milano 
32  Piazza  Leonardo  da  Vinci,  20133  Milano,  Italy 


ABSTRACT 

A  comprehensive  nonsteady  combustion  model,  capable  of  computing 
transient  burning  rates  and  predicting  intrinsic  combustion  stability  of  solid 
rocket  propellants,  is  described  based  on  quasi— steady  thermal  flames.  This 
overall  model  includes  several  appropriate  submodels.  The  transient  flame 
submodel  is  capable  of  describing  spatially  thick  gas— phase  flames,  by 
accommodating  any  flame  thickness  and  allowing  pressure  and/or 
temperature  dependence  of  several  of  the  relevant  parameters.  The  transient 
condensed— phase  submodel  accounts  for  temperature  dependent  properties 
and  possible  volumetric  processes  (distributed  chemical  reactions  and/or 
in— depth  absorption  of  radiation  emitted  by  an  external  source).  The 
pyrolysis  submodel  pictures  irreversible  processes  concentrated  at  the 
burning  surface  and/or  volu  metrically  distributed.  Within  current 
limitations,  the  resulting  overall  transient  combustion  model  is  widely 
applicable,  but  fundamental  work  in  the  area  of  high  temperature  chemical 
kinetics,  characteristic  gas-phase  times,  condensed-phase  heterogeneity,  and 
surface  processes  is  badly  needed. 

The  overall  model  can  either  be  solved  numerically  to  follow  burning 
histories  or  analyzed  on  the  basis  of  stability  theories  to  predict  final  burning 
outcomes.  This  allows  a  convenient  synergism  between  analytical  and 
numerical  methods,  shedding  much  needed  light  if  properly  carried  out. 
Intrinsic  combustion  stability  properties  are  analyzed  from  both  a  linearized 
and  a  nonlinear  but  approximate  viewpoints.  In  the  former  case,  attention  is 
focused  on  sinusoidal  frequency  response  function.  In  the  latter  case,  both 
static  and  dynamic  (or  transitional)  stability  features  are  discussed  by 
introducing  a  nonlinear  algebraic  static  restoring  function.  Dynamic- 
extinction  vs  static  extinction  is  clarified;  deflagration  limits  are  deduced  for 
pressure  and  other  parameters  as  well;  effects  of  combustion  diabaticity  are 
examined;  connections  with  the  Zeldovich— Novozhilov  approach  are  briefly 
discussed.  Bifurcation  plots  are  obtained  suggesting,  for  the  steady— state 
solution,  a  possible  transition  from  time-invariant  to  oscillatory 
self-sustained  for  most  of  the  combustion  parameters;  these  analytical 
expectations  are  numerically  vc’rified.  All  results,  whether  obtained  under 
linear  or  nonlinear  conditions,  may  markc'diy  depend  on  the  enforced 
transient  flame  siilnnodc'l. 


Burning  Behavior  of  Solid  Propellants  at  High  Pressure 
and  High  Loading  Density 


D.  Grune  and  D.  Hensel 

French-German  Research  Institute  of  Saint  Louis  (ISL) 

Saint  Louis,  France 


The  burning  behavior  of  colloid  propellants  is  briefly 
described  by  the  VIEILLE's  law  of  burning,  which  applies  to 
the  range  of  pressure  involved  in  weapons.  For  single-base, 
double-base,  and  heterogeneous  propellants  it  will  be  shown 
by  means  of  closed  vessel  tests  that  the  VIEILLE's  law  is 
not  applicable  to  the  complete  pressure  range.  Therefore, 
the  spatial  law  of  burning  is  modified  by  means  of  an 
empirically  perturbing  term  to  describe  the  combustion  of 
the  foregoing  propellants  as  a  function  of  both  pressure 
and  porosity  of  the  propellant  bed.  To  study  the  parameters 
of  this  novel  spatial  law,  the  pressure  history  in  the 
closed  vessel  test  is  simulated  with  the  interior  ballistic 
code  AMI  -  a  two  dimensional  two-phase  model  -  by  varying 
the  different  parameters,  until  an  approximate  agreement 
between  numerical  and  experimental  investigation  is  rea¬ 
ched  . 
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The  capability  of  numerical  simulations  to  describe  combustion  processes  and  systems 
vdll  always  be  limited  by 

-  Computer  speed  and  memory  available.  For  simuljiting  combustion,  we  need  to 
resolve  extremely  disparate  time  and  space  scales.  Further,  we  will  always  want  to 
resolve  more  and  more  scales  and  contributing  phenomena,  and  therefore  better 
describe  the  system. 

-  Efficiency  of  our  algoritlims  and  their  implementation.  Algorithms  that  make  the 
best  use  of  computer  resourses  usually  imply  some  trade-ofF  between  accuracy  and 
efficiency.  The  important  point  is  to  choose  algorithms  that  are  accurate  enough 
and  most  efficient. 

-  Availability  of  adequate  input  data.  This  refers  to  information  sudi  as  participating 
chemical  species,  chemical  reaction  rates,  thermodynamic  properties,  and  radiation 
intensities. 

This  presentation  addresses  the  first  two  issues,  computers  and  algorithms,  and  con¬ 
siders  the  limits  of  what  we  are  able  to  sim’ilate.  We  have  learned  that  the  way  to 
stretch  this  limit  is  to  break  the  full  problem  into  a  number  of  submodels,  each  repre¬ 
senting  a  contributing  physical  process  .such  as  convection  or  conduction,  solve  for  the 
effects  of  each  process  with  an  accurate  enough  algorithm  optimized  for  the  particular 
computer,  and  then  combine  these  results.  Recently,  we  have  had  access  to  massively 
parallel  computers  in  which  many  computer  processors  carry  out  simultaneous  mathe¬ 
matical  operations  and  then  this  information  is  coordinated.  Such  computers  present 
the  possibility  of  performing  cominitations  {in  order  of  magnitude  or  more  f{vster  than 
conventional  comjiuter.s.  Here  we  describe  numerical  simuhitions  of  several  types  of 
combustion  ])roblems,  ranging  from  laminar  and  diffusion  flames  to  turlmlent  reacting 
flows,  and  use  these  {is  vehicles  for  explaining  the  tnulo-offs  in  accuracy,  efficiency,  and 
inililementation  on  i^arallcl  computers. 
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Dipartimento  di  Ingegnerta  Chtmica,  Universttd  dt  Napolt 
Fiazzale  Tecchio,  80125  NAPOLI,  Italy 

An  unsteady,  two-dimensional  mathematical  model,  including  finite  rate  com¬ 
bustion  kinetics,  buoyancy  effects  and  transport  phenomena  for  both  the  gas  and 
the  liquid  phase,  has  been  numerically  solved  to  simulate  flame  spread  over  liquid 
fuels.  The  vaporization  process  is  described  through  empirical  expressions  relat¬ 
ing  the  vapor  pressure  to  surface  temperature  while  the  the  effects  due  to  surface 
tension  variation  are  properly  accounted  for  by  means  of  the  boundary  conditions 
for  momentum  equations  at  the  liquid /gas  interface. 

In  agreement  with  experiments,  the  flame  propagates  steadily  even  for  initial 
liquid  temperatures  below  the  flash  point.  Moreover,  as  the  initial  liquid  tem¬ 
perature  is  further  decreased,  the  propagation  occurs  according  to  a  pulsating 
behaviour  and,  for  very  low  temperatures,  with  a  uniform  rate  again.  Numerical 
simulations  have  shown  that  oscillatory  phenomena  of  flame  spread  are  always 
associated  with  the  presence  of  surface  tension  driven  motions  in  the  liquid  phase. 
The  presence  of  the  flame  above  a  liquid,  below  the  flash  point,  causes  a  non- 
uniform  temperature  distribution.  Given  the  dependence  of  surface  tension  on 
temperature,  a  gradient  in  the  surface  tension  is  also  established,  which  drives 
motions  from  the  warm  region  under  the  flame  to  the  cool  region  ahead  of  the 
flame.  Thus  an  eddy  is  formed  in  the  liquid  phase,  ahead  of  the  flame,  which 
represents  an  important  mechanism  for  the  pre-heating  of  the  liquid  to  the  flash 
point.  Under  normal  gravity  conditions,  liquid  phase  buoyancy  does  not  seem  to 
affect  the  flame  spread  process.  On  the  contrary,  gas  phase  buoyancy,  because  of 
the  higher  temperatures  due  to  combustion  processes,  plays  a  primary  role  and 
induces  a  convective  flow  field  which  opposes  the  spread  of  the  flame.  As  a  con¬ 
sequence  of  surface  tension  driven  flows  and  of  the  the  shear  at  the  liquid/gas 
interface,  a  circulating  eddy  is  also  formed  in  the  gas  phase  ahead  of  the  flame. 

Gas  and  liquid  phase  convective  heat  transfer  makes  the  surface  temperature 
to  increase  in  a  narrow  region  ahead  of  the  flame,  up  to  give  a  mixture  within  the 
flammability  limits.  The  concurrent  gas  phase  velocities  at  the  flame  leading  edge 
and  the  pre-mixed  flammable  region  over  the  liquid  surface  allows  the  flame  to 
accelerate.  Once  this  high  surface  temperature  region  has  been  passed  across,  the 
flame  decelerates  because  the  liquid,  immediately  ahead  of  the  flame,  is  again  below 
the  flash  point.  Oscillations  in  the  flame  spread  rate  are  associated  with  oscillating 
values  of  the  intensity  of  the  surface  tension  driven  flow.  Surface  velocities  increase 
with  surface  temperature  gradienf.s  which  are  alternatively  large,  when  no  pre- 
mixed  region  exists  ahead  of  th<'  flame,  and  small,  when  the  liquid  ahead  of  the 
flame  has  been  pre-healed  to  the  fla.sh  point.  Tiierefore,  the  formation  and  the 
con.'^umpt ion  of  the  jtre-mixed  (lammable  nsixture  ahea<i  of  the  flame  determine 
the  conditions  for  j)ulsating  sprea<l. 
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Common  approach  to  the  numerical  simulation  of  the 
supersonic  continuous  wave  (CW)  chemical  lasers  is  based  on  successive 
computation  of  the  nozzle  array  and  laser  cavity  flows,  i.e.  on  the  use 
of  the  flow  parameter  values  at  the  nozzles  exit  plane,  obtained  by  an 
independent  calculation  of  fuel  and  oxidant  nozzle  flows,  as  the  boundary 
conditions  for  the  cavity  flow  computation.  Such  approach  apparently 
does  not  permit  to  describe  the  upwind  influence  of  the  cavity  processes 
on  the  nozzle  flow,  which  in  reality  may  be  of  great  significance.  In 
particular,  without  taking  this  influence  into  consideration  it  is 
impossible  to  describe  adequately  the  processes  in  the  recirculation  zone 
behind  the  nozzle  array  base,  which  affect  strongly  the  evolution  of 
mixing  layers  and  combustion  in  the  cavity,  and  consequently  control  the 
laser  efficiency.  Moreover,  the  regimes  of  laser  operating  may  exist, 
when  the  great  difference  of  the  nozzle  exit  pressures  causes  the 
boundary  layer  separation  inside  the  nozzle  array.  In  such  a  situation 
the  combustion  may  start  directly  in  the  nozzles  because  of  diffusion  of 
reactants  into  the  separation  region.  The  only  way  one  could  describe  all 
these  phenomena  is  the  coupled  Navier-Stokes  simulation  of  the  flow  in 
the  domain  of  complex  geometry,  which  contains,  along  with  the  cavity 
region,  at  least  the  supersonic  parts  of  the  fuel  and  oxidant  nozzles. 

In  the  present  work  for  this  purpose  a  new  Navier-Stokes  code  is 
created,  which  is  based  on  recently  developed  Compressibility  Scaling 
Method  (CSM).  The  most  important  feature  of  this  method  is  insensitivity 
of  its  convergence  rate  to  the  flow  Mach  number  Mo,  including  arbitrary 
small  values  of  the  parameter.  Due  to  this  fact,  when  computing  the  flows 
under  consideration  (mixed  subsonic  and  supersonic  flows  with  large  zones 
of  recirculation  and/or  low  Mach  numbers),  CSM  has  significant  advantages 
over  traditional  implicit  methods  for  the  compressible  Navier-Stokes 
equations,  which  efficiency  decreases  sharply  at  low  Mo  values.  When 
simulating  the  combustion  and  nonequilibrium  processes  in  the  HF  chemical 
laser  cavities  CSM  is  coupled  with  time-saving  (without  matrix 
inversions)  numerical  algorithm  of  calculation  of  multicomponent 
diffusive  mass  fluxes  and  with  an  efficient  fully  implicit  method  of 
computation  of  the  radiation  generated  by  the  laser.  Capabilities  of 
proposed  approach  are  illustrated  by  examples  of  its  application  for  the 
■.imulation  of  the  flowficlds,  combustion,  vibrational  relaxation  and 
radiative  transfer  in  HF  CW  chemical  lasers  with  three  different  designs 
of  nozzle  array.  Kinetic  model  of  the  hydrogen-fluorine  combustion,  used 
in  computations,  included  totally  165  elementary  reactions  between  13 
species  (ii,  112,  F,  F2,  He,  DF,  HF(v),  v=0-6  -  vibrational  level  of  the  HF 
molecule).  The  results  presented  demonstrate  both  the  efficiency  of  the 
code  (it  provides  3-5  times  saving  of  \he  computer  time  over  traditional 
algorithms)  and  its  ability  to  predict  complicated  (with  shock  waves, 
recirculating  zones  and  mixing/reacting  layers)  flow  pattern.  The  last 
claim  is  supported  by  good  correspondence  of  numerical  and  experimental 
data  on  flame  structure  in  laser  cavity  and  integral  laser  outpv>t 
characteristics  which  cannot  be  achieved,  oven  qualitatively,  by  the 
common  approach  (successive  computation  of  the  nozzle  and  cavity  flows). 
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A  popular  semi-implicit  scheme,  widely  used  as  an  alternative  to  the  more  expensive  fully- 
implicit  formulation,  is  considered.  The  scheme  is  roughly  equivalent  to  a  pure  explicit 
scheme  in  terms  of  computation  requirements,  is  unconditionally  stable  and  has  the  de¬ 
sirable  feature  of  yielding  non-negative  values  for  species  concentrations.  However,  the 
scliemc  introduces  an  error  resulting  in  chemically-inconsistent  results,  that  are  not  eas¬ 
ily  recognizable  as  such.  These  and  more  properties  of  the  scheme  are  illustrated  and 
discussed,  and  ways  of  overcoming  problems  are  indicated  in  the  present  communication. 
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ABSTRACT 

The  problems  connected  vvitli  the  limitation  of  energy  resources, 
with  the  diversification  of  fuels  and  with  increasing  concern  about  the 
effect  of  emissions  on  the  environment  have  drawn  combustion 
engineers’  attention  to  the  importance  of  improving  combustion 
equipment  design.  At  the  same  lime,  with  computational  abilities 
increasing  each  year,  the  capability  of  obtaining  detailed  information 
about  the  combustion  process,  which  involves  complex  interactions 
between  the  fluid  dynamics,  heat  transfer  and  chemistry,  is  becoming  a 
reality. 

Reviewring  the  available  literature  (Richter,  1985,  Robinson,  1985, 
Troe,  1988,  Bilger,  1980)  wc  realize  that  the  simultaneous  solution  of 
fluid  flow,  heat  transfer  and  complex  kinetics  is  not  yei  sufficiently 
advanced  to  be  applied  to  real  boiler  geometries.  The  aim  of  the 
present  research  is  lo  try  to  develop  a  model  capable  of  coupling  the 
fully  3-D  solution  of  fluid  flow,  heat  transfer  and  complex  chemical 
kinetics  (regardless  of  the  kinetic  scheme  adopted).  In  this  paper 
attention  will  be  focused  mainly  on  the  new  algorithm  for  the  solution 
of  the  kinetic  equations  and  their  coupling  with  the  conservation 
equation  of  the  time-averaged  species. 


Tlie  combustion  inodel  adopted  is  based  on  a  "quasi-global"  scheme, 
made  up  of  a  fuel  demoliiion  reaction  to  CO  and  H2  and  of  a  kinetic 
scheme  for  transformation  to  CO2  and  H2O.  The  assuinpilon  of  "fast 
kinetics”,  generally  utilized  in  the  literature  for  heat  release  calculations 
(e.g.  Boyd  ,1985,  Carvalho  1987)  cannot  be  used  in  tliis  case,  because  it 
permits  us  only  to  calculate  the  final  products  of  the  combustion 
reaction  (CO2,  H2O).  and  not  the  intermediate  species,  that  can 
strongly  influence  NO  formation. 

For  fuel  oil  the  demolition  reaction  to  CO  and  H2  is  assumed  to 
have  intinile  velocity.  For  methane  it  is  considered  with  finite  velocity 
Kp,  as  suggested  by  Malloggi  (1990).  The  kinetic  scheme  for  the 
CO/H2  sy.slem  was  written  considering  the  most  important  reactions 
and  taking  H,  O,  OH  as  the  radical  species.  Twenty-five  reactions  are 
utilized  in  the  present  case,  but  the  program  is  written  in  a  general 
form,  therefore  the  kinetic  scheme  can  be  easily  changed,  allowing  the 
testing  of  different  hypotheses. 

The  species  transport  equations,  after  a  preliminary  transformation 
are  linearized  giving  a  system  of  N  equations  in  N  unknowns  that  are 
solved  for  any  cell  and  for  any  time  step.  This  algorithm  is  very  suitable 
for  the  application  to  a  vectorial  machine,  and  the  solution  implicitly 
satisfies  ma.s.s  conservation,  element  conservation  and  stoichiometric 
relations  with  the  maximum  precision  obtainable  with  the  hardss'are 
utilized,  as  characteristics  of  solutions  in 

As  an  example  of  application  of  the  code,  results  are  presented  of 
the  simulation  of  heat  transfer,  combustion  and  pollutant  formation  in  a 
30  MWt  tangentially-fired  boiler  simulation  facility  where  very  detailed 
in  furnace  measurements  of  temperature  and  main  chemical  species 
were  taken  both  in  standard  ami  staged  combustion  conditions.  Others 
applications  of  the  code  are  de.«;cribcd.  like  design  and  optimization  of  u 
complex  combustion  system  for  ibe  application  of  the  roburning  process 
to  a  160  MWe  power  station  boiler. 
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The  erosion  of  gun  barrels  depends  on  several  influences.  As  we  know  temper¬ 
ature  is  one  of  the  important  quantities.  Erosion  starts  at  the  surface  contacting 
the  hot  propellant  gas.  Till  now  there  exists  an  uncertainty  on  the  real  surface 
temperature  or,  in  other  words,  do  the  experiments  give  the  real  surface  temper¬ 
ature. 

The  approach  of  the  heat  transfer  problem  reported  here  tries  to  stay  in  the 
physics  of  the  event  as  much  as  possible.  The  heat  transport  from  the  hot  pro¬ 
pellant  gas  to  the  gun  barrel  is  computed  by  solving  the  full  Navier-Stokes  equa¬ 
tions  for  the  axisymmetric  in-bore  flow.  This  method  couples  the  center-core  flow, 
the  boundary  layer,  the  heat  transfer,  and  the  heat  conduction  in  the  tube  wall. 
Since  the  boundary  layer  is  considered  the  steep  temperature  gradients  very  close 
to  the  bore  surface  are  available.  Therefore,  the  heat  flux  into  the  wall  is  known 
precisely.  The  numerical  solution  is  based  on  an  ADI  scheme  similar  to  that  of 
Briley-McDonald.  A  one-equation  turbulence  model  completes  the  set  of  equations. 
Comparison  with  temperature  measurement  shows  good  agreement.  The  compu¬ 
tational  method  also  allows  to  include  coated  wall  surfaces. 
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The  up-to-date  nonlinear  high  order  TVD  shock  capturing 
schemes  provide  an  efficient  tool  for  the  simulation  of  the  ideal 
gas  flows  with  complicated  shock-wave  structure.  They  permit 
not  only  to  ensure  the  high  accuracy  of  the  numerical  solution  in 
the  regions  of  it  smooth  behavior  but  also  to  resolve  the  shock 
waves  and  other  modes  of  discontinuities  without  flow  property 
spurious  oscillations.  These  features  of  the  TVD  schemes  can  be 
hardly  overestimated  and  make  them  very  attractive  for  the 
numerical  simulation  of  the  gas-phase  detonation.  However  some 
non-trivial  problems  emerge  when  extending  the  TVD  schemes, 
proposed  for  the  calculation  of  the  flows  of  thermally  and 
calorically  perfect  gases,  to  the  general  case  of  nonequilibrium 
reactive  flows  of  the  real  gas  mixtures.  One  of  such  problems  is 
connected  with  the  development  of  the  appropriate  methods  for 
averaging  of  the  flux  Jacobian  matrices,  which  are  needed  for  the 
approximate  Riemann  solvers  being  an  important  element  of  the  TVD 
schemes  based  on  the  flux-difference-splitting  approach.  The 
methods  of  averaging  proposed  by  now  are  rather  time  consuming  and 
besides  utilize  some  a  priori  non-obvious  assumptions,  which  may 
lead  to  decreasing  of  the  accuracy  of  the  numerical  solution.  In 
the  present  work  a  new  method  of  averaging  is  developed  and 
embedded  into  well  known  Roe's  approximate  Riemann  solver.  It  is 
quite  simple,  provides  exact  satisfying  of  the  U-property, 
formulated  by  P.Roe,  can  be  easily  incorporated  into  any  TVD 
scheme,  and  applicable  for  the  general  case  of  multicomponent 
nonequilibrium  (finite-rate  chemically  reacting)  mixtures  of 
thermally  perfect  gases  with  arbitrary  temperature  dependence  of 
species  internal  energies  and  specific  heat  capacities.  The  above 
properties  of  the  proposed  approach,  coupled  with  high  order  TVD 
scheme  by  S.R.Chakravarthy,  are  illustrated  by  the  two  examples  of 
its  application  for  the  numerical  simulation  of  detonation  waves. 
The  first  example  deals  with  the  oblique  standing  detonation  wave 
holding  by  the  wedge  in  the  supersonic  flow  of  the  stoichiometric 
hydrogen-oxygen  mixture  diluted  by  argon,  and  the  second  one  - 
with  the  evolution  of  the  detonation  induced  by  the  local 
impulsive  energy  addition  in  the  enclosure  filled  with  the  same 
mixture.  In  both  cases  the  same  kinetic  model  of  the  finite  rate 
hydrogen-oxygen  combustion  is  used,  which  includes  9  elementary 
reactions  between  6  species  (h,  0,  OH,  HjO,  H, ,  O,)  . 
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Nuiiiorical  simulation  of  high  speed  reactive  flows  must  face  tlie  cliallenges  prompted  by  the  large 
ditroroiu'i^  witich  may  arise  between  tlie  fluid  dynamic  time  scales,  ry,  and  those  of  the  nonequilibrium 
|)li(>nomena,  r,..  A  nie;isurc  of  the  relative  magnitude  of  these  scales  is  provided  by  the  Damkbhier  number 
=  Ty/r,-),  which  has  a  direct  itifluence  on  the  thickness  of  the  nonequilibrium  layere.  In  the  limiting 
regimes  of  nearly  frozen  {Da  — >  0)  or  nearly  equilibrium  flow  {Da  — >  oo),  there  may  exist  regions  in 
the  flow  domain,  such  as  solid  boundaries  or  shock  waves,  where  very  thin  relaxation  or  chemical  layers 
develo]>.  An  ins\ifficient  numerical  resolution  of  these  layers,  by  reducing  the  accuracy  and  the  stability 
limits  of  the  calculation,  may  affect  the  w'hole  flowfield  by  altering  the  propagation  speed  of  shocks 
loghoter  with  the  ,strurture  of  the  nonequilibrium  layers  at  the  .solid  boundaries. 

With  the  proposed  numerical  model  the  correct  propagation  shock  speed  in  reacting  gases  is  achieved 
by  explicitly  tracking  the  shocks  and  then  by  integrating  the  field  equations  in  quasilinear  form  [1,2],  The 
stability  properties  of  the  scheme  are  enhanced  by  a  semi-implicit  treatment  of  the  stiff  source  terms. 
The  proper  level  of  resolution  at  the  nonequilibrium  layers  is  provided  by  means  of  an  Adaptive  Mesh 
Holinement  (AMR)  technique  (2j. 

The  motivations  backing  the  choice  of  a  tracking  technique  are  briefly  summarized  in  the  following. 

The  great  majority  of  papers  dealing  with  the  numerical  simulation  of  compressible  reacting  flows 
have  adopted  conservative  schemes  developed  in  the  context  of  inert  gasdynamics.  However,  there  seems 
to  exist  some  evidence  (4,5.G!  that  wlien  in  the  nonhomogeneous  system  of  equations  the  reactive  source 
terms  hocome  stiff,  one  of  the  most  appealing  properties  of  tiie  weak  numerical  solutions  may  be  lost; 
the  ability  of  imoviding  the  correct  propagation  speeds  of  the  shocks.  The  origin  of  this  failure  is  related 
to  the  unavoidable  pre,sonco  of  unpbysical  states,  introduced  by  the  shock  capturing  methods  in  order 
to  (lesi  vihe  the  jiropagation  of  shocks,  which  can  ononoously  trigger  the  stiff  source  terms.  In  order  to 
'■deactivate”  the.se  slates  s.ime  kind  of  tracking  seems  to  become  necessary  for  shock  capturing  methods 
as  well. 

Thus,  we  have  deemed  lu'oper  to  follow  a  cUtssical  “fitting”  proceduie  [7j  in  order  to  iisses.s  its  nieriLs 
•  iml  limits  in  the  context  of  reaetive  flows.  Among  the  merits  it  is  immediately  pos.sil)le  to  enlist  the 
following.  The  enforcement  of  jumi)  relations  across  the  filled  shocks,  if  interpreted  as  a  subseale  model 
lilteriiig  out  tlie  relaxation  of  chemical  processes  which  evolve  on  the  shortest  scales  inside  the  sliock 
transit  ion,  can  ])rovide  an  accurate  to(d  to  further  reduce  the  overall  degret'  of  stiffness  dis|)layed  hy  the 
flow  I'oi  I'xaiiqile.  the  "litling”  offers  a  consistent  way  of  computing  the  Jump  relations  according  to  the 
‘'liarliallv  dispersed”  shock  model  [>s]  for  which  the  chemistry  is  frozen  across  the  shock,  and  vibration 
i.'-.  .i.s  the  tbermodynamii  conditions  would  suggest,  alternatively  either  frozen  or  in  equilibrium. 

On  llic  other  baud,  the  issue  of  aii  ellicieul  acliievment  of  .suitable  re.soluliotis  at  the  uonequilibrium 
l.i\eis  li.i-  been  pursued  by  means  of  an  .AMR  tecluiitjue  pre.senled  ill  j2i .  This  teclinique  belongs  to  the 
cl.i.ss  of  embedded,  struct ured  grid.s  anil  it  has  been  designeil  especially  to  obtain  time  accurate  solutions. 

The  generation  of  the  embedded  siibgriiis  is  performed,  indei)en<lent ly  of  the  relineinent  criteria,  once 
.Old  for  all  al  the  beginning  of  the  calculation  according  to  a  completely  recursive  ulgorilhin.  I'he  tnx>  of 
Mibgri'l--  recursivelv  gciieiated  takes  its  root  in  a  relatively  coarst',  structured,  uniform,  orthogonal  grid 
ceiiiaming  a  iiuinbei-.  A',.,//,  of  cells  wbicli  is  a  multiple  of  an  integer  power  of  two.  This  root  grid  i.s 
[i.nlilMneil  inlo  an  e\en  mmibci  of  e<)Ually  sized  grids,  Af;„.,r,  ))reset  by  the  user.  When  any  rclinemelit 
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<’i  itoriuiii  (lolccUs  till!  nci;<l  of  an  iiicroiLscd  rosohitioii  for  a  coll  beloiigiiif;  to  one  of  the  jjarts  of  the  root 
Riiii,  all  cells  of  this  part  have  their  size  halved,  originating  in  this  way  a  new  snbgrid  with  twice  tlie 
resolution  of  the  root  gi  id.  This  procedure  holds  for  all  the  suhgrids  of  the  tree,  'riiercforc,  once  the  pair 
(Nrrti,N,,art)  is  givoii,  tlic  goomctry  and  the  size  of  all  subgrid.s  is  unicjncly  determined,  independently 
of  the  refinement  criterinm  and  of  the  solution  evolution.  The  refinement  critcrium  “simply”  acts  as  a 
mean  to  activate  or  deactivate  the  integration  over  the  already  generated  subgrids. 

'I'ho  time  di.scrotization  and  the  order  of  integration  arc  defined  following  the  asymptotic  idea  of 
kcojiing  the  solution  frozen  on  the  coarser  subgrids  while  letting  it  evolve  on  the  finer  ones.  The  use  of 
structured  grids,  together  with  increasing  the  grid  resolution  by  consecutive  halvings  only,  reduces  to  the 
minimum  the  number  of  interpolations  required  to  let  the  grids  follow  the  projiagating  signals,  and  allow 
an  easier  coordination  of  the  truncation  errors  among  grids  of  different  resolutions.  Both  effects  improve 
the  space  time  accuracy  of  the  calculation. 

The  space-time  link  between  adjacent  subgrids  is  realized  by  adding  extra  computational  points  inside 
the  (hyperbolic)  domain  of  dependence  at  both  sides  of  each  subgrid.  These  extra  points  allow  tlie 
definition  of  a  purely  initial  value  problem  for  each  subgrid,  which  can  thus  be  solved  “independently” 
from  the  integration  on  the  other  subgrids.  This  circumstance  could  become  useful  in  the  context  of 
parallel  computing. 

Once  the  accuracy  requirements  have  been  satisfied  by  the  AMR  strategy,  then  the  degree  of  gasdy- 
namic/chcmical  stiffness  in  terms  of  local  Damkohler  number  is  actually  reduced  or  even  eliminated.  At 
this  point,  without  problems  of  stiffness,  a  fully  coupled  approach  for  the  solution  of  the  system  of  equa¬ 
tions  may  become  feasible  and  perhaps  even  more  advantageous  than  others  based  on  operator  splitting. 
In  fact,  the  fully  coupled  nonlinear  system  can  be  solved  very  efficiently  taking  full  advantage  of  the  high 
convergence  rate  of  Newton-type  methods. 

To  summarize,  in  this  work  the  numerical  integration  of  reactive  flows  is  carried  out,  over  an  adaptive 
computational  domain,  on  the  basis  of  a  fully-coupled,  quasi-linear  formulation  of  the  unsteady,  reactive 
Euler  equations.  The  numerical  method  presented  is  obtained  by  generalizing  the  second-order  A  scheme 
[7]  with  explicit  “fitting”  of  shocks  to  the  case  of  reactive  flows. 

Results  will  be  presented,  in  the  ID  case,  for  //2/O2  combustion,  A  detailed  chemical  kinetics  mech¬ 
anism  is  adopted.  Preliminary  tests  [3],  conducted  with  more  or  less  tight  refinement  criteria,  show  tliat 
the  use  of  the  proposed  AMR  technique,  in  the  context  hereabovc  specified,  requires  a  computational 
w'ork  which  can  be  (>-10%  smaller  than  that  needed  by  a  uniform  grid  of  the  same  resolution.  Similar 
savings  can  be  obtained  in  tenns  of  CPU  time.  Tables  (1)  and  (2)  refers  to  the  following  data. 

Data  for  a  Uniform  Grid 
Number  of  intervals  ;  2018 
Typical  number  of  iterations  :  'lOOO 

Typical  work  (number  of  operations)  :  2  x  2048  x  4000  =  10,384,000  (for  a  two  level  scheme) 

Typical  CPU  time  :  48  hours  of  IBM  RISC  6000/550 

Data  for  the  AMR 
Base  grid  resolution  :  128  intervals 
Number  of  partitions  ;  4 

Finest  level  of  resolution  available  :  equivalent  to  a  uniform  grid  of  2048  intervals 
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Table  1.  Total  work 

Test  Case 

AMR  work  %  Work 

Medium  band  criterium 
Wide  band  criterium 

885,421 

1,526,694  9.3% 

Table  2.  Total  CPU  time 

Test  Case 

AMR  CPU  time  %  CPU 

Medium  band  criterium 
Wide  band  criterium 

180^  6.2% 

309’  10.7% 

Stability  of  travelling  waves  in  deflagration  to  detonation 
transition  models 


Alessandra  Lunardi 

Dipartimento  di  Mateniatica,  Universita  di  Cagliari 
Via  Ospedale  72,  09124  Cagliari,  Italy 


VVe  present  here  some  results  about  the  behavior  near  travelling  wave 
solutions  of  a  class  of  free  boundary  parabolic  problems  with  interface  condi¬ 
tions.  Such  problems  arise  as  specific  models  in  one-dimensional  combustion 
phenomena. 

Our  method  can  be  applied  to  a  number  of  equations  and  systems,  such 
as  the  deflagration-detonation  transition  model,  introduced  by  J.  Ludford 
and  his  collaborators; 

Ui  =  Ur,  +  j;  e  R,  t?  /  ^{t) 

,  *^(4(0,0  =  K]^=f(()  =  -1  /js 

u(-oo,<)  =  0,  u(-foo,f)  =  tioo,  ' 

ti(x,0)  =■  uq(x),  ijEK, 

lYavelling  waves  are  solutions  such  that  ^  is  linear,  say  ^{t)  =  d,  tt 
is  of  the  form  «(<,»?)  =  (/o(»?  +  ct).  After  several  transformations,we  reduce 
problem  (1)  to  an  evolution  equation  in  a  suitable  Banach  'pace,  in  such 
a  way  that  orbital  stability/instability  of  the  TW  solution  is  equivalent  to 
stability/instability  of  the  null  solution  of  the  abstract  equation.  Due  to  the 
interface  conditions,  the  final  equation  is  fully  nonlinear,  i.e.  of  the  form 

r  u'(f)  =  Lt/(<)-f  F(«(0),  <>0  ps 

I  «(0)  =  Uo, 

where  L  i.s  a  linear  operator  in  a  Banach  space  A',  and  the  nonlinear  smooth 
perturbation  F  is  defined  on  the  domain  D{L)  of  L  (more  precisely,  on  a 
neighborhood  of  0  in  D(L))  and  has  values  in  yY.  Geometric  theory  of  equa- 
tion.s  such  as  (2)  includes  results  of  local  existence,  uniqueness,  regularity, 
stability,  liopf  bifurcation  In  particular,  the  principle  of  linearized  stability 
holds,  so  that  the  problem  of  determining  the  stability  of  the  TW  solution 
of  problem  ( 1 )  is  reduced  to  the  problem  of  finding  out  the  spectrum  of 
/  In  our  case.  L  is  nothing  but  the  realization  of  a  second  order  elliptic 
ditb'rential  operator  with  suitable  boundary  conditions  in  the  space  A,  and 
it  is  possilde  to  characterize  its  spectrum.  As  a  consetpience,  we  are  able  to 
giv<'  a  rigorous  mathematical  justification  to  previous  formal  computations. 


Critical  Conditions  for  Detonation  initiation 
by  Nonuniform  Hot  Pockets  of  Reactive  Gases 


H.E.  Longting  and  P.  Clavin 
Laboratoire  de  Recherche  en  Combustion,  URA  1117 
CNRS  &  Universite  d ' Aix-Marseille  I 
Marseille,  France 


A  numerical  study  of  the  critical  conditions  for  detonation 
initiations  in  a  cold  reactive  gaseous  mixture  is  presented 
for  one  dimensional  geometry.  This  study  concerns  initia¬ 
tion  by  a  hot  reactive  pocket  with  nonuniform  initial 
temperature.  The  Zeldovich  criterion  (1980)  for  the  initia¬ 
tion  of  detonation  at  high  temperature  is  found  to  be  well 
verified.  For  a  given  initial  temperature  gradient,  a  CJ 
detonation  is  spontaneously  developed  at  the  temperature 
predicted  by  this  criterion  which  is  typically  about  1100 
K  for  a  characteristic  length  scale  of  few  centimeters.  But 
a  new  phenomenon  is  exhibited:  the  same  initial  temperature 
gradient  yields  a  quenching  of  the  detonation  at  a  lower 
critical  temperature.  Thus,  the  Zeldovich  criterion  is 
shown  to  not  be  sufficient  to  determine  the  critical  condi¬ 
tions  for  detonation  initiation  in  cold  mixtures.  These 
critical  conditions  are  determined  here  in  the  planar 
geometry  for  an  initial  pocket  (slab)  of  hot  reactive 
mixture  with  a  linear  gradient  of  initial  temperature 
embedded  in  a  cold  uniform  mixture.  For  ordinary  reactive 
mixtures,  the  typical  critical  size  of  such  initial  pockets 
to  initiate  detonation  in  cold  gases  at  room  temperature  is 
of  the  order  of  5  cm.  For  a  smaller  size  the  detonation  is 
quenched  before  having  reached  the  room  temperature  and  the 
detonation  initiation  occurs  only  for  larger  sizes.  An 
explanation  of  this  quenching  mechanism  is  provided  as  well 
as  an  analytical  expression  of  the  initiation  criterion. 


Numerical  Simulation  in  Combustion: 


Industry  Needs  and  Perspectives 


A.  Trovati 
FIAT  Avio 
Torino,  Italy 


Even  if  combustion  is  a  very  important  subject  for  gas 
turbine  manufacturers,  this  topic  is  often  still  approached 
quite  empirically  and  the  use  of  numerical  codes  is  rela¬ 
tively  limited. 

Nevertheless,  a  wider  utilization  of  computer  simulations 
would  be  very  effective  in  improving  the  design  process, 
allowing  to  shorten  the  testing  phase,  to  evaluate  a  larger 
number  of  alternative  design  options  and,  eventually,  to 
save  time  and  money. 

The  present  industrial  approach  to  combustion  problems, 
with  particular  reference  to  FIAT  Avio  experience,  is 
described  and  manufacturers'  needs  are  highlighted,  in 
order  to  give  a  guide-line  for  future  cooperation  between 
researchers  and  industries. 

This  cooperation  should  results  in  a  mutual  benefit, 
providing  a  fecund  exchange  of  data  and  experiences.  Tests 
on  technological  rigs  and  data  coming  from  engine  opera¬ 
tions  may  be  a  very  useful  tool  to  check  and  improve  the 
numerical  models  while,  on  the  other  hand,  reliable  models 
would  give  a  dramatic  enhancement  of  design  activities,  as 
already  happened  in  other  fields,  like,  for  example,  stress 
analysis  or  aerodynamics. 


Combustion  St2ibility  of  Solid  Propellants 
ZN-  and  FM-  Approaches 


Boris  V.  Novozhilov 
Institute  of  Chemical  Physics 
Russian  Academy  of  Sciences 
Moscow,  117977,  Russia 


Theoretical  investigations  of  nonsteady  burning  of  solid 
propellants  are  based  mainly  on  the  approximation  which 
takes  into  account  thermal  relaxation  of  the  solid  phase 
only.  Within  this  approximation,  there  are  two  approaches 
to  analysis  named  as  Zeldovich-Novozhilov  (ZN)  and  the 
Flame  Modeling  (FM)  approaches.  The  main  aim  of  this  report 
is  to  point  out  at  the  equivalency  of  the  two  methods.  As 
an  example,  the  problem  of  combustion  stability  is  taken. 
The  well-known  flame  models  such  as  Denison-Baum  (DB)  and 
Krier-T' ien-Sirignano-Summerfield  (KTSS)  models  are  discus¬ 
sed  and  compared  with  ZN-approach.  It  is  shown  that  the 
stability  condition  are  determined  by  only  two  parameters, 
which  characterize  the  dependencies  of  burning  rate  and 
surface  temperature  on  initial  temperature  in  a  steady- 
state  regime. 
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stability  and  Transient  Radiation-Driven 
Solid  Propellant  Combustion 


V.E.  Zarko  and  A.B.  Kiskin 
ICKC,  Novosibirsk  630090,  Russia 


Stability  of  radiation-driven  steady-state  combustion  has 
been  studied  by  the  method  of  small  perturbations.  Stabili¬ 
ty  boundaries  are  determined  depending  on  the  value  of 
radiant  flux  and  propellant  transparency.  The  stability 
domain  is  found  to  become  wider  with  increasing  transparen¬ 
cy  and  radiant  flux. 

Linear  analysis  of  a  burnrate  response  to  harmonically 
oscillating  radial  flux  predicts  a  decrease  of  resonance 
frequency  when  propellant  transparency  and  radiant  flux 
increase.  For  a  resonance-type  burnrate  response^  the 
maximum  phase  shift  is  realized  at  a  frecpiency  less  than 
the  resonance  one.  At  low  frequencies,  the  phase  shift 
value  is  positive  (shift  advance)  for  semitransparent 
propellants  and  negative  for  highly  transparent  ones.  At 
resonance,  the  phase  shift  is  negative  with  absolute  value 
being  less  than  Tf /2 .  Numerical  calculations  revealed  a 
good  agreement  with  the  results  of  linear  analysis. 

Nonlinear  oscillating  burnrate  behavior  is  studied  numeri¬ 
cally.  It  is  shown  that,  due  to  asymmetric  form  of  the 
burnrate  response,  there  is  a  misfit  of  phase  shift  deter¬ 
mined  in  estrema  (max,  min)  of  the  burnrate.  The  effect 
vanishes  with  increasing  propellant  transparency. 

Transient  burnrate  behavior  in  the  case  of  stepwise  radiant 
flux  action  depends  on  propellant  transparency  and  close¬ 
ness  of  a  burning  regime  to  the  linear  stability  boundary. 
If  the  final  state  (self-sustaining  combustion)  is  close  to 
the  boundary  of  stability,  the  burnrate  behavior  after 
radiant  flux  cut  off  is  characterized  by  an  elongated 
transition  for  transparent  propellant  and  extinction  for 
semitransparent  one.  The  reason  is  the  nonconformity 
between  the  enthalpy  excess  in  radiation  driven  combustion 
wave  and  that  in  steady-state  self-sustaining  combustion 
wave . 


Acoustic  Instedsilities  of  Flcuaes 


P.  Clavin 

Laooratoire  de  Recherche  en  Combustion,  URA  1117 
CNRS  &  Universite  d' Aix-Marseille  I 
Marseille,  France 


The  excitation  of  acoustic  vibrations  by  heat  is  related  to 
a  general  phenomenon:  the  thermoacoustic  instability  which 
was  first  pointed  out  by  Rijke  in  the  middle  of  the  XIXth 
century.  Because  of  gas  expansion,  heat  release  in  a  gas 
flow  acts  as  a  volume  source  which  may  easily  excite 
acoustic  vibrations  of  the  cavity.  Lord  Rayleigh  gave  a 
general  criterion  for  a  spontaneous  amplification  of 
acoustic  modes  by  any  local  heat  source:  the  local-energy 
addition  or  substraction  must  be  positively  correlated  with 
the  pressure  fluctuations.  But  this  does  not  tell  us  when 
and  how  such  a  phenomenon  may  be  produced  in  combustion  and 
what  are  the  basic  mechanisms  responsible  for  this  coheren¬ 
ce.  Different  examples  are  presented  in  this  talk.  Recent 
results  of  analytical  studies  concerning  planar  flames  and 
cellular  flames  in  premixed  gas  will  be  first  recalled.  Two 
other  analytical  results  will  be  presented.  A  first  one 
describes  a  new  response  mechanism  of  combustion  of  homoge¬ 
neous  solid  propellants.  It  concerns  the  high  frequency 
range  where  the  non  steady  effects  of  the  gas  phase  become 
dominant.  A  last  one  concerns  the  case  of  flames  propaga¬ 
ting  in  sprays  or  particle-laden  gases  where  inertia  and 
drag  of  particles  may  lead  to  a  strong  acoustic  combustion 
instability.  These  analytical  results  have  been  obtained  in 
the  framework  of  the  asymptotic  analysis  of  non  steady 
flame  structures  in  the  singular  limit  of  an  infinitely 
large  reduced  activation  energy.  Comparison  with  existing 
experimental  data  will  also  be  presented. 


Studies  on  Mathematical  Models  of  Solid  Propellant  Combustion  at  NPV 

Jinqi  Yin 

Northwestern  Polytechnical  University 
Xi'an,  People's  Republic  of  China 


Three  aspects  of  studies  on  mathematical  models  of  solid 
propellant  combustion  conducted  at  Northwestern  Polytechnical 
University  (NPU)  are  reviewed  in  this  paper. 

First,  based  on  SEM  examination  and  X-ray  photoelectron 
spectroscopy  (XPS)  test  on  quenched  surface  of  AP-PU  solid 
propellant,  a  mathematical  model  of  propellant  combustion  is 
developed.  The  model  considers  the  effect  of  molten  binder  covering 
the  AP  crystal  surface,  condensed  phase  reaction,  and  opposed 
gasification  of  AP  covered  with  molten  binder.  In  the  model,  the 
burning  surface  of  solid  propellant  are  divided  into  two  regions,  in 
region  I,  the  AP  crystal  surface  is  covered  with  molten  binder,  and 
the  opposed  gasification  of  AP  is  proposed  to  take  place  under  the 
covering  binder  layer.  In  region  II,  the  BDP  model  well  known  is 
used.  The  two  regions  are  connected  with  the  same  final  adiabatic 
flame  temperature.  The  ratio  of  area  of  region  I  to  that  of  region 
II  and  its  variation  with  pressure  increasing  are  determined 
experimentally  with  XPS  test.  The  model  not  only  can  be  used  to 
explain  the  plateau,  mesa  and  normal  burning  features,  but  also  can 
be  used  to  analyze  the  effects  of  initial  temperature  and  AP 
particle  size  on  burning  behaviors  of  composite  solid  propellants. 
Moreover,  the  model  should  be  capable  of  providing  a  basis  for 
studying  the  evosive  burning  and  unsteady  combustion  of  composite 
solid  propellant,  including  those  with  negative  pressure  exponents. 

Second,  a  simolified  theoretical  model  is  developed  to  consider 
the  ef^fect  of  thermal  radiative  heat  transfer  of  combustion  products 
on  the  combustion  of  aluminized  propellants.  From  the  model,  a 
semiempirical  method  is  deduced  to  predict  the  correlation  of 
propellant  burning  rates  in  a  combustion  bomb,  a  subscale  motor  and 
a  fullscale  motor.  The  method  possesses  the  advantage  of  being 
simple  in  calculation  and  of  requiring  less  experiments. 

The  third  is  about  the  theoretical  model  of  liquid  quench  of  solid 
rocket  motor.  In  the  model,  the  coupling  effects  of  motor 
ballistics,  transient  burning  of  solid  propellants,  liquid  drop 
evaporation  and  heat  transfer  during  cooling  by  liquid  jet  impinged 
on  the  burning  surface  are  considered  in  detail.  The  spacewise  thick 
flame  model  is  used  to  describe  the  transient  burning  behavior  of 
solid  propellants.  With  the  application  of  the  model,  the  liquid 
quench  process  of  solid  rocket  motor  is  predicted  theoretically.  The 
predicted  critical  injection  pressure  drop  of  liquid  jet  and  liquid 
quantity  required  for  reliable  extinction  are  in  good  agreement  with 
those  measured  experimentally. 

In  addition  to  the  aspects  above,  the  theoretical  studies  on 
ignition  and  erosive  burning  of  solid  propellants  are  also  ongoing 
projects  at  NPU. 


Mathematical  Modeling  of  Unsteady  Combustion 
of  Energetic  Solids  Subject  to  External  Radiation 


Quinn  Brewster  and  Steven  Son 
Department  of  Mechanical  and  Industrial  Engineering 
University  of  Illinois  at  Urbana-Champaign ,  USA 


A  theoretical  investigation  is  being  conducted  on  the 
unsteady  combustion  of  energetic  materials  and  solid 
propellants  under  dynamic  radiant  heating  conditions.  The 
analysis  has  been  carried  out  under  the  assumptions  of 
quasi-steady  gas  and  surface  reaction,  homogeneous  propel¬ 
lant,  and  one-dimensional  flame  (or  QSHOD  theory;  "t^ 
approximation"  in  Russian  literature) .  The  classical  linear 
QSHOD  theory  has  been  generalized  to  include  the  effect  of 
external  radiant  flux  with  in-depth  absorption.  This  has 
been  done  in  both  the  framework  of  the  Zeldovich-Novozhilov 
(2N)  phenomenological  theory  and  the  flame  modeling  (FM) 
approach.  The  relationship  of  these  two  approaches  has  been 
demonstrated  for  the  generalized  radiant  flux  case.  The 
nonlinear  QSHOD  theory  has  also  been  extended  to  include 
external  radiation  with  in-depth  absorption  and  some 
numerical  studies  of  nonlinear  effects  have  been  conducted. 
These  studies  have  illustrated  an  effect  of  the  mean 
radiant  flux  on  the  combustion  response  (both  linear  and 
nonlinear)  which  has  previously  been  neglected  by  most 
researchers  in  this  area.  The  theoretical  results  have  also 
illustrated  the  importance  of  in-depth  absorption  for  most 
wavelength-material  combinations.  The  effect  of  external 
radiation,  including  in-depth  absorption,  on  the  intrinsic 
stability  of  steady  deflagration  has  also  been  investigat¬ 
ed.  It  has  been  shown  that,  similar  to  the  case  with 
surface  absorption,  external  radiation  with  in-depth 
absorption  may  have  either  a  stabilizing  or  de-stabil izing 
effect,  depending  on  how  the  temperature  sensitivity  varies 
with  temperature.  The  principle  of  equivalence  of  external 
radiation  and  initial  temperature  has  also  been  re-exa¬ 
mined.  Finally,  analysis  is  extended  to  include  effects 
such  as  distributed  reaction  in  the  condensed  phase  and  a 
fully  unsteady  gas  phase;  some  preliminary  results  are 
presented . 
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Froqiiency  rcs^)onso  functions  of  solid  rocket  propellants,  burning  under 
sinnsoidal  fluctuations  of  ijressurc  or  external  racliant  flux,  arc  discussed  in 
the  framework  of  linear  perturbations.  The  7  flame  model  developc’d  by  this 
research  group  for  spacewise  thick  flames  is  taken  as  reference  for  obtaining 
the  needed  jDcrturlx'd  boundary  condition  at  the  burning  surbice  from  the 
gas— pha.se  side.  Ail  gas-phase  cievelopments  are  expressed  in  terms  of  the 
rharacteristir  gas— phase  time  parameter  projK^sed  by  MTS  modelers  in  1969. 
This  formulation  allows  to  recover  as  special  cases  all  the  results  so  far 
available  in  the  competent  technical  literature  and,  in  addition,  is  a  good 
starting  point  for  future  developments  permitting  the  effects  of  temperature 
on  the  gas-phase  processes  to  lx?  taken  into  account.  Moreover,  this  unifying 
treatment  is  immediately  derivable  from  the  corresponding  nonlinear  analysis 
Ix'ing  carried  out  by  this  res<?arch  group. 

The  alternative  ZN  (Zeldovich— Novozhilov)  method,  mainly  developed 
in  Russia,  is  quickly  touched  upon  and  compared  to  the  FM  (Flame 
Modeling)  method  pursued  in  this  instance.  Uncler  proper  assumptions,  tlie 
two  methods  yield  rather  close  results.  Yet,  considering  the  basic  difficulties 
still  facing  modelers  in  this  technical  arena,  it  is  suggested  that  Ixith  methods 
should  lx?  developed,  hand-in-hand,  in  order  to  .syner^stically  learn  from 
Ixitli. 


'File  (*ffe<  ts  of  average  radiant  flux  intensity  and  arbitrary  surface  layer 
absorption,  previously  neglectcxl  by  this  research  group,  are  included  and 
compared  with  results  recently  published  in  the  open  literature,  d'he  effects 
of  pressure  and  several  other  operating  parameters  are  examined  as  well. 
Among  physical  properties,  the  heat  release  at  the  burning  surface  and 
activation  energy'  of  the  surface  pyrolysis  are  most  influential.  The  intrinsic 
combustion  stability  limits  for  linear  burning  arc  deduced  and  compared  with 
.some  noiilineiu'  expectations.  Finally,  the  problem  of  sharp  names  and 
relevance  of  pressure  de[x?ndent  surface  pyrolysis  is  di.sciissed. 

Fri'sent  di'velopments  are  focu.sed  on  the  influence  of  gas— nha.si? 
temperature  on  fluctuating  heat  feedback  for  both  linear  and  nonlinear 
frecjiu'iiry  r<‘sponse  functions. 


Derivation  of  Simplified  Models  in  Chemical  Kinetics 

by 

D.  A.  Qoussis  S.  H.  Lam 

Mechanical  Engineering  Department  Mechanical  and  Aerospace  Engineering  Department 

Patras  University  Princeton  Univetsity 

Rio-Patras,  65000  Princeton,  NJ  08544 

Greece  USA 


When  massively  complex  chemical  kinetics  problems  are  numerically  investigated,  the 
conventional  methods  of  deriving  simplified  models  arc  no  longer  viable.  In  general,  these 
methods  consist  of  two  steps.  First,  the  intermediates  and  fast  elementary  reactions  must  be 
identified.  Then,  the  steady-state  and  partial-equilibrium  assumptions  must  be  applied  to  the 
intermediates  and  fast  elementary  reactions,  respectively. 

This  process  has  many  shortcomings.  First  of  all,  the  difficulty  in  identifying  the 
intermediates  and  fast  reactions  is  proportional  to  the  size  of  the  full  mechanism.  Furthermwe, 
the  success  of  the  conventional  analysis  requires  a  large  amount  of  computations  and  relies 
heavily  on  the  intuition  and  experience  of  the  investigator.  Finally,  the  results  of  such  an 
analysis,  depend  on  the  assumptions  made  for  the  initial  and  operating  conditions. 

The  Computational  Singular  Perturbation  (CSP)  meth^  is  a  general  method  for  the 
treatment  of  nonlinear  ODE'S  which  exhibit  a  boundary-layer  type  of  stiffness.  As  a  result,  the 
method  is  very  useful  in  the  analysis  of  chemical  kinetics  problems.  CSP  identifies  the  fast 
reactions  and  intermediates  and  derives  the  appropriate  relations  among  the  elementary  reaction 
rates,  i.e.  the  steady-state  and  partial-equilibrium  approximations.  Since  CSP  is  a 
computational  method,  the  size  and  complexity  of  the  kinetics  mechanism  as  well  as  the 
dependence  of  the  result.s  to  the  initial  and  operating  conditions  pose  no  limitations  on  the 
successful  use  of  the  method. 

The  relations  among  different  elementary  reaction  rates  impose  some  constrains  on  the 
evolution  of  temperature  and  the  species  mass  fraction.  What  really  controls  the  system's 
evolution  within  these  constrains  is  not  clear.  To  answer  this  question,  the  original  system  of 
differential  equations  modelling  the  physical  process  must  be  simplified  by  the  appro]xiate  use 
of  the  relations  among  different  elementary  reaction  rates.  As  it  will  be  shown  later  on,  this  is 
an  enormous  task  for  an  investigator,  especially  if  the  kinetics  mechanism  is  large. 

CSP  can  numerically  perform  this  task  requiring  no  outside  input  of  any  kind.  The  ability  of 
CSP  to  properly  identify  the  elementary  reaction  rates  which  locally  control  the  evolution  of 
temperature  and  species  mass  fractions  and  to  simplify  the  system  of  differential  equations 
provides  important  informations  in  the  analysis  of  chemical  kinetics  mechanisms.  First,  the 
minimum  set  of  unknowns  can  be  selected  so  that  no  effort  is  wasted  in  monitoring 
unimportant  to  the  investigator  dependent  parameters.  In  addition,  the  proper  kinetics  constants 
can  be  adjusted  so  that  the  numerical  model  provides  a  good  representation  of  the  experimental 
data. 

The  power  of  the  CSP  method  in  acquiring  all  the  appropriate  physical  insights  and  in 
simplifying  a  large  kinetics  mechanism  will  be  demonstrated  for  the  case  of  homogeneous 
methanol  oxidation  at  constant  pressure. 


International  Workshop. 


on 

MATHEMATICAL  METHODS  IN  COMBUSTION 


List  of  participants 


Centro  di  Culture  Scientifica  "A.  Volta" 
Villa  Olmo,  Como,  Italy 


18-22  May  1992 


UST  OF  WORKSHOP  PARTICIPANTS  WITH  COMPLETE  ADDRESS 
Villa  Olmo,  Como.  Italy,  1S-20  May  1992 


Mr.  Salvatore  BARUFFO 
IRC/CNR 

328  Via  Diodeziano 
80124  Napoli,  Italy 
Telex  : 

Phone  :  (0039-81)  570-5864 
Fax  :  (0039-81)  762-2915 

Prof.  M.  Quinn  BREWSTER 

Department  of  Mechanical  and  Industrial  Engineering 

University  of  Illinois  at  Urbana-Champaign 

144  Mechanical  Engineering  Building 

1206  West  Green  Street 

Urbana,  IL  61081,  USA 

Telex; 

Phone;  (001-217)  244-1628 
Fax  ;  (001-217)  244-6534 

Dr.  Vladimiro  CALARESE 

EOARD,  USAF 

223  Old  Marylebone  Road 

London  NW1  5TH,  UK 

Telex; 

Phone:  (0044-71)  409-4423 
Fax  :  (0044-71)  724-1433 

Prof.  Giancarlo  CHIATTI 

Dipartimento  di  Meccanica  ed  Aeronautica 

Universitii  di  Roma  "La  Sapienza* 

18  Via  Eudossiana 
00184  Roma,  Italy 
Telex: 

Phone:  (0039-6)  4458-1 
Fax  :  (0039-2)  488-1759 

Prof.  Paul  CLAVIN 

LRC  Laboratoire  de  Recherche  en  Combustion 

URAD  1117  CNRS,  Universit6  de  ProverKe 

Centre  Saint  Jerome 

13397  Marseille  Cedex,  France 

Telex;  402876  FACSTJE 

Phone:  (0033-9)  128-8124 

Fax  :  (0033-9)  163-5261 

Mr.  Giovanni  COLOMBO 
Dipartimento  di  Energetics 
Pditecnico  di  Milano 
32  Piazza  Leorrardo  da  Vinci 
20133  Milano,  Italy 
Telex:  333467  POUMI  I 
Phone:  (0039-2)  2399-3871 
Fax  :  (0039-2)  2399-3940 


Dott.  Tarcisio  COLOMBO 
Alfa  Romeo 
FIAT  Auto 

20020  Arese,  Ml.  Italy 
Telex  : 

Phone  ;  (0039-2)  9339-2972 
Fax  ;  (0039-2)  938-0769 

Dr.  Gaetano  CONTINILLO 
IRC/CNR 

328  Via  Oiodeziano 

80124  Napoli,  Italy 
Telex:  722392  INGENA  I 
Phone:  (0039-811  570-5864 
Fax  :  (0039-81)  762-2915 

Mr.  Fabio  COZZI 
Dipartimento  di  Energetica 
Politecnico  di  Milarx) 

32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex:  333467  POUMI  I 
Phone:  (0039-2)  2399-3872 
Fax  :  (0039-2)  2399-3940 

Dr.sa  Roberta  DAL  PASSO 
lAC/CNR 

Istituto  ApplicazionI  Calcolo  'Mauro  Picone' 

Consiglio  Nazionale  delle  Ricerche 

137  Viale  del  Polidinico 

00161  Roma,  Italy 

Telex: 

Phone:  (0039-6)  8847-0262 
Fax  :  (0039-6)  440-4306 

Dr.  Giuliano  D’ANDREA 
Science  Adviser 
USASETAF 
36100  Vicenza.  Italy 
Telex  : 

Phone  :  (0039-444)  505-409 
Fax  :  (0030^144)  505-409 

Prof.  Luigi  DE  LUCA 
Dipartimento  di  Energetica 
Politecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex:  333467  POUMI  I 
Phone:  (0039-2)  2399-3912 
Fax  :  (0039-2)  2399-3940 

Prof.  Colomba  DI  BLASI 
Dipartimento  di  Ingegneria  Chimica 
University  di  Napoli 
80  Piazzale  V  Tecchio 

80125  Napoli 
Telex: 

Phone:  (0039-81)  768-2264 
Fax  .  (0030-81)  593-6936 


Prof.  Luciano  GALFETTI 
Dipartimento  di  Energetica 
Pditecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex:  333467  POUMI  I 
Phone;  (0039-2)  2399-3926 
Fax  .  (0039-2)  2399-3940 

Mr.  Dietrich  GRUNE 

ISL  (Institut  Franco-Allemand  de  Recherches  de  Saint  Louis) 

Group  B6/PCH,  Dept.  B 

12  rue  de  I’lndustiie,  Boite  postale  No.  34 

F-68301  Saint  Louis,  France 

Telex:  8811386 

Phone:  (0033)  8969-5000  ext.  5073 
Fax  ;  (0033)  8969-5002 

Dr.  Rudi  HEISER 

Fraunhofer-lnstitut  fur  Kurzzeitdynamik 
Emt-Mach-lnstitut  Abteilung  fur  Ballistik 
18  Hauptstrasse 

D-7S58  Weil  am  Rhein,  Germany 
Telex: 

Phone:  (0049-7621)  7-1067 
Fax  :  (0049-7621)  7-6925 

Prof.  Nickolai  KIDIN 

Institute  Chemical  Kinetics  and  Combustion 

USSR  Academy  of  Sciences 

Siberian  Branch 

630090  Novosibirsk,  Russia 

Telex;  133148  KING  SU 

Phone:  (007-3832)  35-7030 

Fax  :  (007-3832)  35-2350 

Prof.  Alessandra  LUNARDI 
Dipartimento  di  Matematica 
Facdt^  di  Scienze 
Universit^i  di  Cagliari 
Via  deirOspedale 
Cagliari,  Italy 

c/o:  Scuda  Normale  Superiore 
Pisa,  Italy 
Telex: 

Phone:  (0039-50)  597-290 
Fax  :  (0039-50) 

Ing.  Fabrizio  MAGUGLIANI 
Convex  Computer  S.p.A. 

Centro  Direzionale  Cdleoni 
Palazzo  Orione,  3 
1 7  Viale  Cdleoni 
20041  Agrate  Brianza,  Ml,  Italy 
Telex:  326423  SIAVBC  I 
Phone:  (0039-39)  605-8024 
Fax  :  (0039-39)  609-1079 


Dr.  David  M.  MANN 

Associate  Director 

US  Army  Research  Office 

Engineering  arxJ  Environmental  Science  Division 

Research  T riangle  pArk 

NC  27709-221 1 ,  USA 

Telex: 

Phone.  (001-919)  549-4249 
Fax  :  (001-919)  549^1348 

Prof.  Bernard  J,  MATKOWSKY 

Department  of  Engineering  Sciences  and  Applied  Mathematics 
McCormick  School  of  Engineering  and  Applied  Science 
Northwestern  University 
Evanston,  Illinois  60208,  USA 
Telex: 

Phone:  (001-708)  491-5396 
Fax  :  (001-708)  491-2178 

Prof.  A.  G.  MERZHANOV 

Institute  of  Structural  Macrokinetics 

Russian  Academy  of  Sciences 

142432  Chemogoloka,  Moscow  Region,  Russia 

Telex:  411701  SFERA  SU 

Phone:  (007-95) 

Fax  :  (007-95) 

Prof.  Boris  V.  NOVOZHILOV 
Institute  of  Chemical  Physics 
Russian  Academy  of  Sciences 
4  Kosygin  Street 
117977  Moscow  V-334,  Russia 
Telex:  411701  SFERA  SU 
Phone:  (007-95)  939-7127 
Fax  :  (007-95)  938-2156 

Dr.  Elaine  S.  ORAN 

Laboratory  for  Computational  Pnysics 

Code  4044/4004 

Naval  Research  Laboratory 

Washington,  DC  20375,  USA 

Telex: 

Phone:  (001-202)  767-2960 
Fax  :  (001-202)  767-4798 

Prof.  Carlo  Domenico  PAGANI 
Dipartimento  di  Matenatica 
Pditecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Mtano,  Italy 
Telex:  333467  POLIMI  I 
Phone:  (0039-2)  2399-4517 
Fax  .  (0039-2)  2399-4568 

Ing.  Sauro  PASINI 
ENEL-CRTN 
Via  A.  Pisano,  1 20 
56100  Pisa,  Italy 
Telex 

Phone:  (0039-50)  53-5630 
Fax  (0039-50)  53-5651 


Dr.ssa  Francesca  PISTELLA 
lAC/CNR 

Istituto  Appticazioni  Calcdo  "Mauro  Picone" 

Consiglio  Nazionale  delle  Ricerche 

137  Viale  del  Pdidinico 

00161  Roma 

Telex; 

Phone;  (0039-6)  8847-01 
Fax  ;  (0039-6)  440-4306 

Anna  Maria  PULLE' 

Dipartimento  di  Energetica 
Pditecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex;  333467  POUMI  I 
Phone;  (0039-2)  2399-3801 
Fax  ;  (0039-2)  2399-3940 

Ing.  Alberto  QUARANTA 

Direttore 

CNPM/CNR 

Istituto  Propulsione  ed  Energetica 
Consiglio  Nazionale  delle  Ricerche 
69  Viale  Francesco  Baracca 
20068  Peschlera  Borromeo,  Ml,  Italy 
Telex; 

Phone;  (0039-2)  715015 
Fax  ;  (0039-2)  714063 

Prof.  Alfio  QUARTERONI 
Dipartimento  di  Matematica 
Pditecnico  di  Milano 
32Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex;  333467  POUMI  I 
Phone;  (0039-2)  2399^525 
Fax  ;  (0039-2)  2399-4588 

Ms.  Letizia  RATTI 
Dipartimento  di  Energetica 
Politecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex;  333467  POUMI  I 
Phone;  (0039-2)  2399-3872 
Fax  ;  (0039-2)  2399-3940 

Dr.  Roy  E.  REICHENBACH 
Aeronautics  and  Mechanics  Branch,  Chief 
European  Research  Office 
USARDSG  (UK) 

223  Old  Marylebone  Road 
London  NW1  5TH,  UK 
Telex: 

Phone:  (0044-71)  409^23 
Fax  ;  (0044-71)  724-1433 


Prof.  Sandro  SALSA 
Dipartimento  di  Matematica 
Politecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex;  333467  POUMI  I 
Phone:  (0039-2)  2399-4553 
Fax  :  (0039-2)  2399-4568 

Mr.  Ugo  SCAGUONE 
Dipartimento  di  Energetica 
Politecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex;  333467  POUMI  I 
Phone;  (0039-2)  2399-3872 
Fax  ;  (0039-2)  2399-3940 

Prof.  Michael  L  SCHUR 

Department  of  Physics  and  Mechanics 

State  Technical  University 

197198  Saint  Petersburg,  Russia 

Telex;  121345  PTB  SU 

Phone:  (007-812)  238-5353 

Fax  :  (007-812)  232-2027 

Prof.  Michael  Kh.  STRELETS 
State  Institute  of  Applied  Chemistry 
1 4  Dobrolyubov  Av. 

197198  Saint  Petersburg,  Russia 
Telex;  121345  PTB  SU 
Phone:  (007-812)  238-5353 
Fax  :  (007-812)  232-5844 

Prof.  Alberto  TESEI 

Direttore 

lAC/CNR 

Istituto  Applicazioni  Calcolo  'Mauro  Picone" 

Consiglio  Nazionale  deile  Ricerche 

137  Viale  del  Polidinico 

00161  Roma 

Telex: 

Phone.  (0039-6)  8847-01 
Fax  :  (0039-6)  440-4306 

Ing.  Alberto  TROVATI 

Fiat  Avio  Direzione  Progettazione 

112  C.  so  Ferrucci 

10138  Torino,  Italy 

Telex: 

Phone;  (0039-11)  3302-802 
Fax  :  (0030-11)  3302-486 

Ing.  Mauro  VALORANI 

Dipartimento  di  Meccanica  e  Aeronautica 

Universit^i  di  Roma  'La  Sapienza" 

18  Via  Eudossiana 

00184  Roma,  Roma,  ITALY 

Telex: 

Phone:  (0039-6)  4458-5281 
Fax  :  (0039-6)  488-1759 


Prof.  Maurizio  VERRI 
Dipartimento  di  Matematica 
Politecnico  di  Milano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex:  333467  POUMI  I 
Phone;  (0039-2)  2399-4526 
Fax  :  (0039-2)  2399-4568 

Ing.  Luigi  VIGEVANO 

Dipartimento  di  Ingegneria  Aerospaziale 

Politecnico  di  Milano 

40  Via  Golgi 

20133  Milano,  Italy 

Telex;  333467  POUMI  I 

Phone;  (0039-2)  2399-4057 

Fax  :  (0030-2)  2399-4034 

Prof.  Jinqi  YIN 

Department  of  Astronautical  Engineering 
Northwestern  Polytechnical  University  (NPU) 
P.O.  Box  164 

Xi'an,  Shaanxi  Province  710072,  China 
Telex.  70185  NWTRUY  CN 
Phone;  (0086-29)  53351 
Fax  ;  (0086-29)  751-959 

Prof.  Anna  ZARETTI 
Dipartimento  di  Matematica 
Politecnico  di  MPano 
32  Piazza  Leonardo  da  Vinci 
20133  Milano,  Italy 
Telex;  333467  POUMI  I 
Phone:  (0039-2)  2399-4512 
Fax  :  (0039-2)  2399^568 

Prof.  V.E.  ZARKO 

Institute  of  Chemical  Kinetics  and  Combustion 

USSR  Academy  of  Sciences 

Siberian  Branch 

630090  Novosibirsk,  USSR 

Telex.  133148  KING  SU 

Phone:  (007-3832)  35-7030 

Fax  :  (007-3832)  35-2350 


